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Abstract We study non-viscous and viscous holographic dark energy models for a homogeneous and 
isotropic flat Friedmann-Robertson-Walker Universe in f(R,T) gravity. We find that the Hubble 
horizon as an 1R cut-off is suitable for both the models to explain the recent accelerated expansion 
of the Universe. The cosmological parameters like deceleration parameter and statefinder parameters 
are discussed in each model. In non-viscous model a constant deceleration parameter is found which 
shows that there is no phase transition. The constraints on the parameters are obtained to analyze the 
fixed point values of statefinder parameters of SCDM and ACDM models. We know that the phase 
transition is required to explain the accelerated expansion of the Universe and this is possible if both 
the parameters would be time-dependent. Therefore, we extend our analysis to viscous holographic 
dark energy model to investigate whether this viscous model with the same IR cut-off could be helpful 
to find the phase transition. We find that this model gives a time-dependent deceleration parameter 
which achieves a smooth phase transition of the Universe. We also find the time-varying statefinder 
pair which matches with ACDM model. We plot the trajectories in r — s and r — q plans to discriminate 
our model with the existing dark energy models and obtain the quintessence like behavior for the 
suitable values of parameters. 
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1 Introduction 

It is strongly believed that the Universe has entered a phase of the accelerated expansion which has 
been confirmed by the recent observations like supernovae la [1, 2], cosmic microwave background 
radiation [3, 4], baryon acoustic oscillation [5] and Planck data [6]. Within the framework of general 
relativity (GR), the cause of the acceleration can be attributed to the existence of a mysterious com¬ 
ponent of the Universe dubbed as “dark energy” (DE), which makes up ~ 70% of the total cosmic 
energy in the Universe. The ACDM model presents the simplest and most successful description of the 
recent accelerated expansion scenario and accommodates the observations very well. Despite of many 
attractive features, it has some theoretical problems like fine-tuning and cosmic coincidence problems 
[7-9]. To overcome from these problems, a number of dynamical dark energy models such as scalar 
field (quintessence, phantom, k-essence, etc.) models [10-13], chaplygin gas models [14], holographic 
dark energy (HDE) models [15-17], etc. have been explored in the literature . 

In the recent years, the HDE models have been emerged as a viable candidates to explain the 
problems of modern cosmology. The HDE models explain the recent accelerated expansion as well as 
the coincidence problem of the Universe [18, 19]. The concept of HDE is based on the holographic 
principle proposed by’t Hooft [20] and found it’s roots in the quantum field theory. Cohan et al. [21] 
have shown that in the quantum field theory, the formation of black hole set a limit which relates 
UV cut-off A to IR cut-off L. According to the authors, the quantum zero-point energy ph = A 4 of 
a system of size L should not exceed the mass of a black hole of the same size, i.e., L 3 ph < LM 2 , 
where M p = (87tG) - 1 / 2 is the reduced Planck mass. In a paper [18], Li has taken the largest allowed 
L to saturate this inequality and thus obtained dark energy density of the Universe pu = 3c 2 M 2 L~ 2 , 
known as HDE density. In the formalism of HDE, the Hubble horizon is a most natural choice for the 
IR cut-off, but it leads to a wrong equation of state (EoS) of dark energy [22]. However, Pavon and 
Zimdhal [19], and Banerjee and Pavon [23] have shown that the viable EoS of dark energy could be 
achieved by taking the interaction between HDE and dark matter (DM). 

On the other hand, the modified theories of gravity such as f(R) gravity [24, 25], /(G) gravity [26, 
27], f(R, G ) gravity [28, 29], etc. have also been proposed to explain the recent accelerated expansion 
of the Universe. The f(R) gravity is one of the simplest and successful modified theories of GR, which 
fits with the observations very well. Recently, Harko et al. [30] have proposed a new modified gravity 
theory known as f{R , T) gravity, where R as usual stands for the Ricci scalar and T denotes the trace 
of energy-momentum tensor. This modified theory presents a maximal coupling between geometry 
and matter. A number of authors [31-36] have discussed the modified f(R,T) gravity in different 
context to explain the early and late time acceleration of the Universe. In a recent paper [37], the 
authors have discussed the viscous cosmology in this theory which shows the recent phase transition 
of the Universe. The HDE models have not been yet discussed in detail in the framework of f(R,T) 
gravity. In some papers [38, 39], reconstruction of f(R, T) gravity from HDE and anisotropic model 
of HDE have been discussed. 

In this paper, we are interested to discuss HDE with Hubble horizon as an IR cut-off in the frame¬ 
work of f{R,T) gravity. It has been shown in GR and Brans-Dicke theory that the Hubble horizon 
as an IR cut-off is a suitable candidate to explain the recent accelerated expansion if we consider 
interaction between HDE and DM [19, 23]. As it is known, the f(R,T) gravity has inbuilt interaction 
between geometry and matter, therefore, it will be interesting to discuss HDE with Hubble horizon as 
an IR cut-off in this modified theory. The observations from different probes suggest that the DM is 
a non-interacting content of the Universe. Therefore, it is assumed that only HDE component of the 
total matter (HDE +DM) interacts with the geometry of the Universe i.e. T of /(i?, T) is the trace 
of the energy-momentum tensor of HDE only. We show that the Hubble horizon as an IR cut-off is 
suitable to explain the recent accelerated expansion in this modified gravity theory. 

The Hubble parameter H and the deceleration parameter q are well known cosmological parame- 
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ters which explain the evolution of the Universe. However, these two parameters can not discriminate 
among various DE models. In this context, Sahni et al. [40, 41] have introduced a new geometrical 
diagnostic pair {r, s}, known as statefmder parameters, which is constructed from the scale factor 
and its derivatives up to the third order. The statefmder pair {r, s} is geometrical in the nature as 
it is constructed from the space-time metric directly. Therefore, the statefmder parameters are more 
universal parameters to study the DE models than any other physical parameters. In a flat ACDM 
model, the statefmder pair has a fixed point value {r, s} = {1,0}. One can plot the trajectories in 
r — s and r — q planes to discriminate various DE models. We discuss the statefmder diagnostic and 
obtain the fixed point values of statefmder pair {r, s} = {1,1} and {r, s} = {1,0} as in the case of 
SCDM and ACDM models, respectively, under suitable constraints. 

To be more realistic, the prefect fluid Universe is just an approximation of the viscous Universe. 
The dissipative processes in the relativistic fluid may be modeled as bulk viscosity. The phenomenon 
of the bulk viscosity arises in the cosmological fluid when the fluid expands (contracts) to fast due 
to which the system is out of thermal equilibrium. Then, the effective pressure become negative to 
restore the thermal equilibrium [42]. Therefore, it is natural to consider the bulk viscosity in an accel¬ 
erating Universe. It has been shown that inflation and recent acceleration can be explained using the 
viscous behavior of the Universe, and plays an important role in the phase transition of the Universe 
[43-50]. The concept of viscous DE has been discussed extensively in the literature [51-53]. Feng and 
Li [54] show that the age problem of the Ricci dark energy can be alleviated using the bulk viscosity. 
Motivated by the above works, we extend our analysis to viscous HDE with the same IR cut-off which 
gives the recent phase transition of the Universe. We obtain the statefmder parameters for viscous 
HDE which achieve the value of ACDM model and show consistency with the quintessence model for 
suitable value of parameters. 

The paper is organized as follows. In the next section we discuss the formalism of f(R, T) gravity 
theory and present its field equations. In section 3 we discuss the non-viscous HDE model and find 
the exact power-law solution of the scale factor which avoids the big bang singularity. We also find 
the cosmological parameters like deceleration parameter and statefmder parameters and discuss their 
behaviors. Section 4 describes the viscous HDE model and its solution. Section 5 presents the sum¬ 
mary of our findings. 

2 The formalism of modified f(R,T) gravity theory 

The general form of the Einstein-Hilbert action for the modified f(R, T) gravity in the unit 8ttG = 1 
is as follows [30, 36]: 

S=\f d*xyf—g [f(R, T) + 2 C m ], (1) 

where g stands for the determinant of the metric tensor g^ u , R is the Ricci scalar and T represents the 
trace of the energy-momentum tensor, i.e., T = T)}', while C m denotes the matter Lagrangian density. 
The speed of light is taken to be unity. As usual the energy- momentum tensor, of matter is 
defined as 

rp _ ^ S(^/~gC m ) . . 

1 j 

In fact, this modified gravity is the generalization of f(R) gravity and is based on the coupling between 
geometry and matter. The corresponding field equations have been derived in metric formalism for 
the various forms of f(R,T). 

Varying the action (1) with respect to the metric tensor for a simple form of f(R, T ) = R+f(T), 
i.e., the usual Einstein-Hilbert term plus an f(T) correction [30, 36] which modifies the general 
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relativity and represents a coupling with geometry of the Universe, we get the following field equations. 

~ 2 RQiiv ~ ~ {Rfiv + ©jui/)/ (T) + — f{T)g liI/ , (3) 

where a prime denotes derivative with respect to the argument. The tensor 0 M1/ in (3) is given by 

d 2 C 

= -2 + g^Cm - 2 3 Q/3 QgpvQgap • ( 4 ) 

The matter Lagrangian C m may be chosen as C m = —p [30], where p is the thermodynamical pressure 
of matter content of the Universe. Now, Eq. (4) gives = —2 T^ v —pg^ v . Using this result, Eq. 
(3) reduce to 

R/j-v — 2 Rd^v = T^v + {Rfiv + P g^v)f (T) + 2 f{ r R)9tivi (5) 

which are the held equations of the modified f(R,T) gravity theory. 

Here, we are interested to study the behavior of HDE in this modified theory for a spatially ho¬ 
mogeneous and isotropic hat Friedmann-Robertson-Walker (FRW) spacetime, which is expressed in 
comoving coordinates by the line element, 

ds 2 = dt 2 — a 2 {t){dx 2 + dy 2 + dz 2 ), (6) 

where a{t) stands for the cosmic scale factor. In what follows we study the non-viscous and viscous 
HDE models with deceleration parameter and statehnder parameters in f(R,T) gravity theory to 
describe the recent acceleration. 

3. Non-viscous holographic dark energy cosmology 

In this model, let us consider the Universe filled with HDE plus pressureless DM (excluding 
baryonic matter), i.e., 

rji rph _i_ rpm ( r 7\ 

1 HV — 1 /iv ' 1 iiv i l'/ 

where T^ v and T™ represent the energy-momentum tensors of HDE and DM, respectively. Many 
authors have described the recent accelerated expansion by assuming the interaction between HDE and 
DM in the different theories of gravity. In this paper, instead of taking the interaction between HDE 
and DM to describe the recent acceleration, we consider that the HDE interacts with the geometry of 
f(R,T) gravity. This is due to the fact that this modified gravity theory has the interaction between 
matter and geometry. Therefore, we consider T = g^T^ as the trace of energy-momentum tensor of 
HDE. The generalized Einstein equations (5) yield 

3H 2 = p m + ph + (p m + Ph + Ph) f' (T) + -f(T), (8) 

2H + 3H 2 = -p h + ±f(T), (9) 

where p m , ph and ph denote the energy density of DM, the energy density of HDE and the pressure 
of HDE, respectively. An overdot denotes the derivative with respect to cosmic time t. As the field 
equations (8) and (9) are highly non-linear, therefore, we assume f(T) = aT [see, ref. 30], where a 
is a coupling parameter. Now, the field equations (8) and (9) reduce as 

3 H 2 = p. m + ph + a(p m + ph + Ph) + — Q; T, 


(10) 
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2H + 3H 2 = - Ph + i a T. (11) 

The equation of state (EoS) and the trace of energy-momentum tensor of HDE are given by ph = WhPh 
and T = ph — 3ph, respectively. Now, from (10) and (11), a combined evolution equation for H can 
be written as 

2 H + (1 + Ql)[(l + Wh)p m + ph] = 0. (12) 

In the literature, various forms of HDE ( the general form is ph = 3?M 2 L~ 2 , where c 2 is a 
dimensionless constant, M p stands for the reduced Planck mass and L denotes the IR cut-off radius) 
have been discussed depending on the choices of IR cut-off such as Hubble horizon, future event 
horizon, apparent horizon, Granda-Oliveros cut-off, etc. In this work, we consider the Hubble horizon 
(L = H -1 ) as an IR cut-off to describe the recent acceleration. The corresponding energy density ph 
is given by 

p h = 3 c 2 H 2 . (13) 

Form (10) and (13), the energy density of DM can be written as 


3 (ac 2 Wh — 3 ac 2 — 2 c 2 + 2) 2 

= -2(T+~o)- H ■ 


(14) 


Using (13) and (14) into (12), we finally get 

3 

H + -(Sac 2 ^ + 2 c 2 Wh — a? + 2) H 2 = 0, 


(15) 


which, on solving gives 

H = _I_ 

Co + | (3 ac 2 Wh + 2c 2 Wh — ac 2 + 2 )t’ 

where Cq is a positive constant of integration. Eq. (16) can be rewritten as 


H = 


_ Ho _ 

1 + ^p(3 ac 2 Wh + 2c 2 Wh — ac 2 + 2 )(t — t 0 ) ’ 


(16) 


(17) 


where Hq is the present value of the Hubble parameter at the cosmic time t = to, the time where the 
HDE starts to dominate. Using the relation H = -. the cosmic scale factor a is given by 


a = Ci 


1 + -Ho(3ac 2 Wh + 2c 2 Wh - ac 2 + 2 )(t - t 0 ) 


3(.3ac 2 w h +2c 2 -w h -ac 2 + 2) 


(18) 


where Ci is an another positive constant of integration. One can rewrite a as follows 


a = a o 


1 + —Hg(3ac 2 Wh + 2 c 2 Wh — at? + 2 )(t — t 0 ) 


3(3oc- i » t + 2c- ! « k -ac- ! + 2) 


(19) 


where a o is the present value of the scale factor at the cosmic time t = to- We obtain the power-law 
of evolution of the Universe which avoids the big-bang singularity. 

The deceleration parameter q, which is defined as q = —ad/a, is a geometric parameter which 
describes the acceleration or deceleration of the Universe depending on the negative or positive value. 
In this case, the deceleration parameter is given by 

1 3 

q = - + -?(3aw h + 2w h - a). 


(20) 
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Here, we obtain a constant deceleration parameter as expected due to the power-law of the evolution. 
As we observe that for a given value of uy,, one can obtain an accelerated expansion for coupling 
parameter a satisfying the constraint a > 3 c\i-zw h ) • For example, one can take a > — | for c 2 = 1 
and u>h = — 1 to observe the accelerated expansion. Thus, the HDE with Hubble horizon as an IR 
cut-off can successfully explain the accelerated expansion in the framework of f(R, T) gravity without 
assuming the interaction between HDE and DM in contrast to the works done in [19, 55]. It is to be 
noted that this model does not show the phase transition as the deceleration parameter is constant. 

In order to get a robust analysis to discriminate among DE models, Sahni et al. [40, 41] have 
introduced a new geometrical diagnostic pair {r, s}, known as statefinder parameters, which is con¬ 
structed from the scale factor and its derivatives up to the third order. The statefinder pair {r, s} is 
geometrical in the nature as it is constructed from the space-time metric directly. The statefinder pair 
{r, s} provides a very comprehensive description of the dynamics of the Universe and consequently 
the nature of the DE. It is defined as 


a r — 1 

of/3’ 3(g- 1/2)' 


( 21 ) 


In this model, we obtain the statefinder parameters r and s as 


9 9 

r = — c 4 (3awh + 2 wr — ot ) 2 + — c 2 (3awh + 2 Wh — a) + 1, 
8 4 


( 22 ) 


s = y (3 aw h + 2wh - a) + 1. (23) 

We observe that the statefinder parameters {r, s} are constant and the values of these parameters 
depend on the coupling parameter a, constant c and EoS parameter Wh of HDE. In the papers [40, 
41], it has been observed that SCDM model and ACDM model have fixed point value of statefinder 
pair (r, s} = {1,1} and {r, s} = (1, 0}, respectively. In our work, it is observed that SCDM model can 
be achieved for a = for any values of c. Thus, for a suitable value of wy, which may be obtained 

by observations, we can find the coupling parameter a for which {r, s} = {1,1} and viceversa. 

In a paper, Li et al. [56] have studied the Planck constraints on HDE and obtained the tightest 
and self-consistent value of constant c from Planck+WP+BAO+HST-l-lensing as c = 0.495 ± 0.039. 
Therefore, let us consider here and thereafter c = 0.5 for further discussion which is lying in this 
observed range. Now, assuming any values of Wh, we can get the value of coupling parameter a, e.g., 
c = 0.5 and Wh = —1, we have a = 1.5, which achieves the fixed point of ACDM model (r, s} = {1,0}. 
Observations show that EoS of HDE is not exactly the same of cosmological constant (w/, = —1), in 
fact it may lie in quintessence region (—1 < Wh < —1/3) or phantom region Wh < —1- The coupling 
parameter a provides us the flexibility to obtain the ACDM model for any values of Wh lying in the 
quintessence/phantom region. At the boundary of quintessence region i.e. Wh = —1/3 and c = 0.5, 
we have a = 11/3 for which (r, s} = {1, 0}. Here, we observe that as we decrease the value of Wh from 
-1/3 to -1, the value of a decreases from 11/3 to 3/2 to obtain (r, s} = {1, 0}. Similarly, one can get 
{r, s} = {1,0} in phantom region for suitable value of a. 


4 Viscous holographic dark energy cosmology 


In non-viscous HDE model we have obtained the constant value of deceleration parameter which 
does not describe the phase transition. But, the astronomical observations show that the phase 
transition is an integral part of the evolution of the Universe. Therefore, in this section, it would 
be of interest to investigate whether a viscous HDE with the Hubble horizon as an IR cut-off could 
be helpful to find the phase transition, i.e., time -dependent deceleration parameter and statefinder 
parameters in order to elucidate the observed accelerated expansion of the Universe. 



7 


In an accelerating Universe, it may be natural to assume that the expansion process is actually a 
collection of state out of thermal equilibrium in a small fraction of time due to the existence of possible 
dissipative mechanisms. In an isotropic and homogeneous FRW model, the dissipative process may 
be treated via the relativistic theory of bulk viscosity proposed by Eckart [57] and later on pursued 
by Landau and Lifshitz [58]. It has been found that only the bulk viscous fluid remains compatible 
with the assumption of large scale homogeneity and isotropy. The other processes, like shear and heat 
conduction, are directional mechanisms and they decay as the Universe expands. Bulk viscosity can 
produce an accelerated expansion even without dark energy matter due to the presence of an effective 
negative pressure. Recently, the present authors [37, 49] have studied the effect of viscous fluid in 
f(R,T ) gravity and discussed the recent accelerated expansion of the Universe. 

Using the Eckart formalism for dissipative fluids [57], we can assume that the effective pressure of 
HDE is a sum of the thermodynamical pressure ( ph ) and the bulk viscous pressure (II), i.e., 


P e ff =Ph + P = Ph- 3C H, 


(24) 


where £ is the positive coefficient of the bulk viscosity. Now, the matter Lagrangian is taken as 
£m — - P eff for which Eq. (4) gives 0 M „ = — 2T fJjV — P e f/g In this model we follow the same 
concept as discussed in Section 3 to analyze the behavior of the Universe. Thus, we assume that the 
viscous HDE matter interacts with the geometry of the Universe. Using f(T) = aT, Eq. (3) yields 
the field equations for viscous HDE in the framework of f(R, T) gravity as 

3 H~ = p m + ph + a(p m + Ph + Ph — 3 (H) + —aT, (25) 

2 H + 3 H 2 = -p h + 3 (H + i aT . (26) 

Using T = ph — 3(ph — 3 (H) into (25) and (26), a single evolution equation of H is given by 

2 H+(a + 1 )(p m + p h + p h ) - 3(a + l)(H = 0. (27) 

From (25), we get 

Pm = m 3 -n H [(«c 2 ^ - 3ac 2 - 2c 2 + 2 )H - a<]. (28) 

2(a + l) 

Now, Using (13) and (28) into (27), we get 

H + ^(3 ac 2 w h + 2c 2 w h - ac 2 + 2 )H 2 - |(3a + 2 )(H = 0. (29) 


Eq. (29) is solvable for H if the coefficient of bulk viscosity ( is known. Many authors have studied 
the cosmological models by assuming the various forms of the bulk viscous coefficient (for review, see 
[59]). Here, we assume the bulk viscous coefficient to be a constant, i.e., ( = (o, which is the simplest 
form of ( = Co + (lP [37, 60] by taking = 0. Now, Eq. (29) reduces to 


H+ 'j{3a c 2 w h + 2 c 2 w h - a c 2 + 2 )H 2 - ^(3 a + 2)C, 0 H = 0. 


The solution of (30) is given by 

H = 


e f(3a+2)Co* 


. (3 a c 2 w h +2 c 2 w h -a c 2 +2) 3f3a+2')<:n* ’ 

L2 ^ (3a+2)Co 


where C 2 is a constant of integration. The scale factor a in the terms of cosmic time t is 

„ 2 „„ , o „ 2 „„ „,„2 


a = c 3 


C 3ac-w h + 2c w h -ac +2) 3 (3Q+2)Coi ' 

2+ (3a + 2)Co 


3(3ac 2 w h +2c 


yc 2 +2) 


(30) 


(31) 


(32) 



where C 3 > 0 is another constant of integration. Eq. (32) can be rewritten as 


a = ao 


, (3 a?w h + 2?w h - a? + 2) tt 3 (3a+2)Co(t _ 
(3a + 2) Co H ° [e 


to) 


- 1 ) 


3(3occ'^w h +2c' 2 z, 


2 + 2 ) 


(33) 


One can observe that the model avoids the big-bang singularity. In this case, the deceleration param¬ 
eter is given by 


Q = 


[(3ac 2 Wh + 2?w w - ac 2 + 2 )H 0 - (3a + 2)£ 0 ] 


ei 


(3a+2)Co(i-to) 


- 1 . 


(34) 


It is observed that the value of q is time-dependent which comes due to the introduction of bulk viscous 
term in HDE. The phase transition of the Universe can be explained using this value of deceleration 
parameter. The deceleration parameter must change it’s sign from positive to negative to explain 
the recent phase transition (deceleration to acceleration) of the Universe. In fact, q must change the 
sign at the time t = to because we have assumed to is the time where the viscous HDE begins to 
dominate. In other words, the Universe must decelerate for t < t-o (matter dominated epoch) and 
accelerate for t > to (HDE dominated epoch). We observe that the Universe shows the transition 
from decelerated to accelerated phase for a = —1(1 + 3w/ g r )_(° c ;~ 77 ~) at cosmic time to- Therefore, 
the value of coupling parameter a can be obtained for a given value of uy,, which may be obtained 
from the observations, or vice-versa to get the transition. Thus, we have shown that the bulk viscous 
HDE with Hubble horizon as an IR cut-off can explain the recent phase transition of the Universe in 
the framework of f(R,T) gravity. 

Next, we discuss the another geometrical parameters, i.e., statefinder parameters. In this case, the 
statefinder parameter r is obtained as 


——[(3 a?w h + 2 c 2 w h - at? + 2)H 0 - (3a + 2)<o] 

16-fzo 

( 3a?Wh + 2 ?Wh — a? + 2 )Hq — (3a + 2)£o (3 ac 2 Wh + 2?Wh — a? — 2)Hq 


(35) 


H 0 e 


§(3a+2)Co(i-to) 


+ 


(3<*+2) Cob - ^o) 


+ 1 . 


The second statefinder parameter s is not given here due to complexity but one can find it by us¬ 
ing the values of q and r from (34) and (35) in s = o(q?i/ 2 ) • Our m °del reproduces the fixed 
point value {r, s} = {1,0} of ACDM model when the parameter a satisfies the condition a = 
—1(1 + .)• For this value of a, the statefinder pair is independent of time and re¬ 

mains fixed throughout the evolution as in ACDM model. Indeed, we have obtained time-dependent 
statefinder pair which means that a general study of the behaviour of this pair is needed. We plot the 
trajectories in r — s and r — q planes for some particular values of parameters a and Wh to discriminate 
our model with existing models of DE. 

In Figs, la, b, the fixed points {r, s} = {1,1} and {r, s} = {1, 0} have been shown as SCDM model 
and ACDM model, respectively. It is obvious from both the figures that for any values of a and Wh, 
the viscous HDE model always approaches to the ACDM model, i.e., {r, s} = {1,0} in late the time 
evolution. However, in early time of the evolution our model can approach in the vicinity of SCDM 
model for some values of a as can be seen in Figs, la, b. It is interesting to note that there exist some 
negative values of a for which viscous HDE model starts from ACDM in early time and approaches 
to the same ACDM model during the late time of evolution. 

In the quiessence model with constant EoS (Qi-model) [40, 41] and the Ricci dark energy (RDE) 
model [61], it has been shown that the trajectories in r — s plane are vertically straight lines. In the 
both models, s is constant throughout the evolution of the Universe, while r increases in RDE model 
and decreases in Qi-model starting from the initial point r = 1. It has also been shown in [40, 41] 
that the trajectories for the quintessence scalar field model {Q 2 — model) where the scalar potential 
V(4>) varies as V{<j>) oc > 1 approach asymptotically to the ACDM model in the late time. 
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Comparing this viscous HDE model with Q\— model and RDE model, we find that our viscous HDE 
model produces the curved trajectories which approach to ACDM model in the late time. Further, 
we observe that our model almost shows the similar trajectories like Q 2 — model for some values of at 
and Wh in r — s plane, e.g., for a = —0.5, Wh = —1 and a = —0.4, Wh = —0.5 as shown in Figs, la, 
b, respectively, show the similar trajectories as Q 2 —model for /3 = 2 [40, 41]. 




( a ) ( b ) 

Figure 1. The trajectories in r — s plane are plotted for Wh = — 1 in left panel (a) and for = —0.5 in right panel 
(b) for different values of a. Here, we have taken c = 0.5, Hq = 1 and f = 0.05. 




q q 

(a) (b) 

Figure 2. The trajectories in r — q are plotted for Wh = — 1 in left panel (a) and for = —0.5 in right panel (b) for 
different values of a. Here, we have taken c = 0.5, Hq = 1 and £ = 0.05. 

Figs. 2a, b plot the trajectories in r — q plane. The SCDM model and SS model (steady-state 
cosmology) have been shown by the fixed points {r, q} = {1,0.5} and {r, q} = {1,-1} , respectively. 
The horizontal line in Figs. 2a, b represents the evolution of the trajectory corresponding to ACDM 
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model starting from the fixed point {r, q} = {1,0.5} of SCDM model and end at the fixed point 
{r, q} = {1,-1} of SS model. It can be observed from both the figures that for any values of a and 
Wh, the viscous HDE model always approaches asymptotically to the SS model, i.e., {r, q} = {1, —1} 
as ACDM model in late time evolution. However, in the early time of evolution the model passes 
through SCDM model for some values of a as shown in Figs. 2a, b. The trajectories corresponding 
to our model start from the higher values of the pair {r, q} but after a certain time span show the 
(52-model like behaviour. Again, comparing our viscous HDE model with Qi-model, RDE model and 
<52—model, we find that viscous HDE model is compatible with Q 2 — model. 


5 Conclusion 

In GR and Brans-Dicke theory, some authors [22, 55] have found that the Hubble horizon is not a 
viable candidate to explain the accelerated expansion of the Universe. However, Pavon and Zimdahl 
[19], and Banerjee and Pavon [23] have shown that the interaction between HDE and DM can de¬ 
scribe the accelerated expansion. Therefore, it is clear that one can observe accelerated expansion 
if the interaction between the different matter contents is considered. In this work, we have studied 
non-viscous and viscous HDE models with Hubble horizon as an IR cut-off in the frame work of mod¬ 
ified f{R,T) gravity. The f(R,T) gravity theory presents a maximal coupling between geometry and 
matter. Therefore, we have explored the consequences of the coupling of matter with the geometry 
of the Universe instead of taking the interaction between HDE and DM as many authors have stud¬ 
ied. However, we have assumed that only HDE of total matter (HDE+DM) couple with geometry. 
We have investigated the possibility whether the Hubble horizon as an IR cut-off could explain an 
accelerated expansion in f(R , T) gravity by assuming the interaction between HDE and geometry. 
We have shown that the non-viscous and viscous HDE models with Hubble horizon as an IR cut-off 
can explain accelerated expansion in the frame work of this modified theory for suitable values of the 
parameters. Further, we have investigated statefinder pair {r, s} to discriminate our non-viscous and 
viscous HDE models with other existing DE models. We summarize the results of these two models 
as follows: 

In non-viscous HDE model, we have found an accelerated expansion under the constraint of pa¬ 
rameters. In this case, we have obtained constant deceleration and statefinder parameters. Due to 
constant g, it is not possible to analyze the phase transition of the Universe. We have found the fixed 
points {r, s} = {1,1} and {r, s} = {1,0} of SCDM and ACDM model, respectively, for suitable choice 
of the parameters. Thus, non-viscous HDE model is consistent with SCDM and ACDM models. 

In viscous HDE model, we have obtained the recent phase transition of the Universe as the de¬ 
celeration parameter comes out to be time-dependent. In this model, the statefinder parameters are 
the function of cosmic time t. These time-dependent parameters are possible due to the inclusion of 
bulk viscous fluid in HDE model which could explain the recent phase transition in a better way. It 
is interesting to note that the viscous HDE model gives the ACDM model fixed point {r, s} = {1,0} 
and remains fixed in ACDM model throughout the evolution for a specific value of a as discussed in 
section 3. The statefinder diagnostic have been discussed through the trajectories of r — s and r — q 
planes as shown in Figs. 1 and 2 to discriminate our model with the existing DE models. In Figs, 
la, b, it has been observed that some of the trajectories pass through the vicinity of SCDM during 
early time but ultimately all approach to ACDM model in the late time. In Figs 2a, b, it can be seen 
that one of the trajectories passes through SCDM model for a suitable value of a in early time but all 
the trajectories approach to SS model in the late time of evolution. It has been noticed that for some 
values of a the trajectories of the viscous HDE model are similar to the trajectories of <52—model [40, 
41]. Therefore, the viscous HDE in the framework of f(R,T) gravity gives more general results in 
comparison to ACDM and Q 2 —model at least at the level of statefinder diagnostic as we are able to 
achieve the behavior of both the models. 
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